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Abstract 



We propose and study a model of dephasing due to an environment of bistable 
fluctuators. We apply our analysis to decoherence of Josephson Qubits, in- 
duced by background charges present in the substrate, which are also re- 
sponsible for the 1// noise. The discrete nature of the environment leads 
to a number of new features which are mostly pronounced for slowly moving 
charges. Far away from the degeneracy this model for the dephasing is solved 
exactly. 
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Any quantum system during its evolution gets entangled with the surrounding environ- 
ment. This effect is known as decoherence Besides the understanding of its role in 
many fundamental questions, decoherence is studied because it will ultimately limit the 
performance of a quantum computer 0. Solid state nanodevices seem the natural arena to 
fulfill the requirements of large scale integrability and flexibility in the design, though, due 
to the presence of many kinds of low energy excitations in the environment, decoherence 
represents a serious limitation. Proposals to implement a quantum computer using supercon- 
ducting nanocircuits are proving to be very promising and several experiments already 
highlighted the quantum properties of these devices 

In superconducting nanocircuits various sources of decoherence are present as fluc- 
tuations originating from the surrounding circuit, quasiparticle tunneling, fluctuating back- 
ground charges (BC) and flux noise. In this Letter we introduce and study a model for 
decoherence due to a discrete environment, which describes what is considered the most 
serious limitation for Josephson qubits in the charge regime, i.e. the decoherence originat- 
ing from fluctuating charged impurities. The common belief is that this effect originates 
from random traps for single electrons in dielectric materials surrounding the island. These 
fluctuations cause the 1// noise directly observed in Single Electron Tunneling devices P,^. 

The system under consideration is a Cooper pair box [Q. Under appropriate conditions 
(charging energy Eq much larger than the Josephson coupling Ej and temperatures ksT <^ 
Ej ) only two charge states are important and the Hamiltonian of the qubit Tiq reads 

^ 5Ec Ej 

where the charge basis { | ) , | 1 ) } is expressed using the Pauli matrices, and the bias 
6Ec = Ecil—CxVx/ e) can be tuned by varying the applied gate voltage Vx- The environment 
is a set of BCs electrostatically coupled to the qubit. The resulting total Hamiltonian is 

n = nQ- K^^vMk +T.'^<i (i) 

Hi = Eciblbi + YllTkiclibi + h.c] + ^ Ekicl^Cki . 

k k 

Here Hi describes an isolated BC: the operators hi (6-) destroy (create) an electron in 
the localized level Ed- This electron may tunnel, with amplitude Tki to a band described 
by the operators Cki, cl^ and the energies Eki- For simplicity we assume that each local- 
ized level is connected to a distinct band. An important scale is the total switching rate 
■ji = 27iAf{eci)\Ti\'^ {Af is the density of states of the electronic band, jT^jp ^ l^iP); which 
characterizes the classical relaxation regime of each BC. Finally the coupling with the qubit 
is such that each BC produces a bistable extra bias fj. The same model for the BCs has been 
used in Ref. and explains experiments on charge trapping in systems of small tunnel 
junctions. 

Our aim is to investigate the effect of the BC environment on the dynamics of the 
qubit. In this Letter we focus on decoherence during the quantum time evolution, due to an 
environment which produces 1// noise. Thus we are concerned with BCs with a distribution 



of switching rates ~ I/7 , in the range [7^,, 7m] fTT 
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We used several techniques as second order perturbation theory in the couphngs Vi, or 
in the BC-band couphngs T^j, Heisenberg equations of motion and a real-time path-integral 
analysis [p!^-p!4|. The picture which emerges is that the decoherence produced by each BC 
depends qualitatively on the ratio = Vi/'ji, so it is convenient to distinguish between 
weakly coupled BCs {gi <^ 1) and strongly coupled BCs (in the other regimes). We stress 
that a strongly coupled BC does not necessarily have large coupling Vi with the qubit. Indeed 
in the physical situation we discuss in this Letter the energy scale associated to the total 
extra bias produced by the set of BCs is much smaller than the level splitting of the qubit, 
\J5E^ + Ej. To summarize our results, we find that as far as decoherence is concerned, 
a single weakly coupled BC behaves as a source of gaussian noise, whose effect is fully 
characterized by the power spectrum of the unperturbed equilibrium fluctuations of the 
extra bias operator Vi h\hi , given by Si{u) = vf/[2 cosh^(/3eci/2)] 7i/ {■yf + cu^). Each weakly 
coupled BC contributes independently to decoherence. Instead our quantum mechanical 
treatment yields that decoherence due to a strongly coupled BC shows pronounced features 
of its discrete character, as saturation effects when ^fj ^ 1, and dependence on initial 
conditions. 

We consider, for the sake of clarity, two special operation points for the qubit: (i) charge 
degeneracy, 6Ec = and (ii) the case where tunneling can be neglected, Ej = 0. For this 
latter case, where pure dephasing occurs without relaxation, we find an exact solution. 

(i) Charge degeneracy {5Ec = 0) - Relaxation (F^) and dephasing (F<^) rates of the qubit 
given by the Golden Rule [|g,|T§ 

rr = l^?"" = l^iEj) (2) 

depend only on the power spectrum ^(ci;) = J2iSi{uj) at u = Ej. This simple result 
would readily allow an estimate of the rates Tr^,/, [|15| from independent measurements of 
S{uj) P,|5|JTB[|. Being second order in Vi, the Golden Rule is appropriate only for weakly 
coupled BCs. We study the general problem by using the Heisenberg equations of motion. 
For the average values of the qubit observables (o"a), (sum over a, P = x,y, z is implicit) we 
obtain {h = 1) 

N 

{a a) = Ej eap{ap) +riap^Vi {h%ap) 

i 

where N is the total number of BCs and e^^ = —e^y = 1, rj^y = —rjy^ = 1 are the only 
nonzero entries. In the rhs averages of new operators which involve also the localized levels 
and the bands are generated. New equations have to be considered and the iteration of 
this procedure yields an infinite chain. A closed set of 3 + 3A^ equations is obtained by 
factorization of high order averages: we ignore the cumulants {h\hih^jhj)c and {h\hih^jhjaa)c 
for i ^ j and insert the relaxation dynamics for the BCs in the approximated terms. This 
method gives accurate results for general values of gi even if Vi/Ej is not very small, as we 
checked by comparing with numerical evaluation of the reduced density matrix of the qubit 
with few BCs. Results are presented in Figs. where the time Fourier transform of {az{t)), 
proportional to the average charge on the island, is shown. We assumed factorized initial 
conditions for the qubit and the BCs. We first consider a set of weakly coupled BCs in the 
range [10~^, 10] Ej which determine 1/f noise in a frequency interval around the operating 
frequency. The coupling strengths Vi have been generated uniformly with approximately 
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zero average and magnitudes chosen in order to yield the amphtude of typical measured 
spectra (extrapolated at GHz frequencies). These BCs are weakly coupled, and 

determine a dephasing rate which reproduces the prediction of the Golden Rule Eq.(^ 
(Fig.0 dotted line V^^/Ej ^ 1.65 x 10~^). Now we add a slower (and strongly coupled 
^70 = ^^0/70 = 8.3) BC, in order to extend the 1/f spectrum to lower frequencies. The added 
EC gives negligible contribution to the power spectrum at Ej so according to Eq.(0) it 
should not modify F^. Instead, as shown in Fig.|I], we find that the strongly coupled BC 
alone determines a dephasing rate comparable with that of the weakly coupled BCs. The 
overall F^ is more than twice the prediction of the Golden Rule. If we further slow down 
the added BC we find that increases toward values ~ 70. This indicates that the effect of 
strongly coupled BCs on decoherence saturates (we discuss later similar results for the case 
of pure dephasing, where this conclusion can be made sharp). 

In this regime we observe also effects related to the initial preparation of the strongly 
coupled BC (see Fig.^. Finally we checked different sets of BCs with the same power 
spectrum. They yield larger decoherence if BCs with g ^ 1 are present. The Golden 
Rule result Eq. (Q) underestimates the effect of these strongly coupled BCs. All the features 
presented above are a direct consequence of the discrete nature of the environment. 

Pure dephasing {Ej = 0) - In the absence of the tunneling term Eq.(|l]) is a model for 
pure dephasing. The charge on the qubit island is conserved, [H, = 0, but superpositions 
of charge states dephase. Hence the off-diagonal elements of the reduced density matrix 
of the qubit in the charge basis, or equivalently the averages of the raising and lowering 
operators a±, decay in time. 

It is possible to show by direct calculation that for a product initial condition, {a±) fac- 
torizes exactly in independent contributions of each BC, ((T_(t)) = ((T_(0)) Yij exp{—i{6Ec — 
Vj/2)t}fj(t). Using a real-time path-integral technique |]TB|], the general form of fj(t) in 
Laplace space is obtained 

f (X) = X + K,,AX)-iv,l25p] 

where = 1 — 2{b'jbj)t=o- The kernels i^ij(A) and K2j{X) are given by formal se- 
ries expressions. An explicit form is obtained in the Non Interacting Blip Approximation 
(NIBA) [|13|,|ll, ifi,,(A) = 7,-, K2,j{X) = -7,V7r[^(l/2 + /5/27r(A - ^e^,-)) - V^(l/2 + /?/27r(A + 



iecj))]- In order to appreciate the validity of the NIBA result, we notice that it is also 
obtained by the Heisenberg equations with the only assumption that the band is in equilib- 
rium. 

We now discuss the results in the physically relevant limit where the BCs have an inco- 
herent dynamics. In this case an analytic form for (cr_(t)) is found 

= exp {-i5E^t} exp {-F(t) + iE{t)} , (4) 

where we defined Uj = ^1 — g'j — 2i gj tgh {f3ecj/2) and Aj = ^ (l + aj — i gj^ . This 
result can be obtained both by approximating the kernel K2j{X) ~ K2j{0), valid in the 
limit Eci ,Vi ,'-fi <^ KbT, and as the exact result of a semiclassical analysis. In this last case 
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the coupling operator J2iViblbi is substituted by a classical stochastic process £{t), sum of 
random telegraph processes, yielding 



-iSEct 



((exp[-z dt'S{t')])) (5) 



where ((. . .)) is the average over the possible realizations of S{t) with given initial conditions 

The form of Eq.(^) elucidates the different role of weakly and strongly coupled BCs 
in the decoherence process. Dephasing due to each BC comes from the sum of two expo- 
nential terms. If gj -C 1 only the first term is important and the corresponding rate is 
~ l/[4 cosh^(/?ecj/2)] v'^hj, the Golden Rule result. If Qj ^ 1 the two terms are of the same 
order, and the decay rate is ~ 7^, the switching rate of the BC. The main effect of strong 
coupling with the qubit is a static energy shift. That is for the slower BCs {gj ^ 1) the 
contribution to the decoherence rate saturates to ~ 7^. At short times t <^ 1/7^- the initial 
conditions, which may take the values = ±1, determine the transient behavior. 

We now apply this result to sets of BCs which produce 1// noise. We stress that while 
saturation of dephasing due to a single BC is physically intuitive, it is not a priori clear 
whether this holds also for the 1// case, where a large number (~ I/7) of slow fiuctuators 
(strongly coupled BCs) is involved. In Fig.|^ we show the results for a sample with a number 
of BCs per decade = 1000 and with Vi distributed with small dispersion around {v) = 
9.2 X IQ^ Hz. Initial conditions = ±1 are distributed according to = Speq, the 

equihbrium value. To illustrate the different role played by BCs with gj <^ 1 and gj 1, we 
consider sets with •jm = lO^^Hz and different 7^- In this case the dephasing is given by BCs 
with 7j > 10'^ Hz ~ {v)/10. The main contribution comes from three decades at frequencies 
around (v). The overall effect of the strongly coupled BCs (7j < (f)/10) is minimal, despite 
of their large number. 

Finally we compare our results with the gaussian approximation. It amounts in estimat- 
ing the average in Eq.(^ by its second cumulant and taking Sp^ = 6peq 



12 (t) = - du S{u) . (6) 

TT JO LU^ 

This formula fails to describe BCs with gj ^ 1. For instance r2(t) at a fixed t, scales with 
the number of decades and does not show saturation. The gaussian approximation should 
become correct if the environment has a very large number of extremely weakly coupled BCs. 
We check this limit by comparing T2{t) with Eq.(^). The power spectrum S{uj) is identical 
for all the curves in Fig.^ but is obtained by sets of charges with different rid and (v) . The 
gaussian behavior is recovered for t 3> l/7m if is large (all the BCs are weakly coupled). 
If in addition we take 6p^ = 6peg in Eq.(|D, T(t) approaches r2(t) also at short times. Hence 
decoherence depends separately on and (v) , whereas in the gaussian approximation only 
the combination n^iv'^), which enters S{uj), matters. In other words the characterization of 
the effect of slow sources of 1// noise requires knowledge of moments of the bias fluctuations 
higher than ^(ti;). 

In conclusion we studied decoherence due to an environment of bistable charges. We 
found that the average coupling between individual BCs and the qubit is an important scale 
of the problem: BCs such that 7^ <C Vi show pronounced features of their discrete dynamics. 
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as saturation and transient behavior. The physical picture we obtain for the decoherence 
effects due to BCs is not sensitive to the details of the model Hamiltonian (|l]), but is mainly 
determined by the discrete character of the environment. Thus this approach can also be 
applied to different physical systems as phase qubits in the presence of flux noise 0. 

We thank G. Castelli, Y. Nakamura, M. Palma, M. Rasetti, A. Shnirman, M. Sassetti, D. 
Vion, U. Weiss and A. Zorin for discussions. This work has been supported by EU (IST-FET 
SQUBIT) and INFM-PRA-SSQI. 
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FIGURES 




I ^ — 

0.999 1.001 1.003 1.005 

(B/Ej 

FIG. 1. The Fourier transform az{uj) for a set of weakly coupled BCs plus a single strongly 
coupled BC (solid line). The separate effect of the coupled EC {go = 8.3, thin line) and of the set 
of weakly coupled BCs (dotted line), is shown for comparison. Inset: corresponding power spectra. 
In all cases the noise level at Ej is fixed to the value S{Ej)/Ej 3.18 x 10~^. 
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FIG. 2. The Fourier transform 0-2(0;) for a set of weakly coupled BGs plus a strongly coupled 
BC (^0/70 = 61.25) prepared in the ground (dotted line) or in the excited state (thick line). Inset: 
corresponding power spectra (the thin line corresponds to the extra BC alone). 
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FIG. 3. Saturation effect of slow BCs for a 1/f spectrum. Relevant parameters 
{{v) = 9.2 X 10^ Hz, rid = 1000) give typical experimental measured noise levels and reproduce 
the observed decay of the echo signal in charge Josephson qubits. Couplings Vi are distributed 
with dispersion {Av)/{v) = 0.2. T{t) is almost unaffected by strongly coupled charges (the label is 
the number of decades included). 
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FIG. 4. Ratio T{t)/r2(t) for a 1/f spectrum between 7^ = 2 x 10^ and 7m = 2 x 10^ with 
different numbers of BCs per decade: (a) = 10^, (b) = 4 x 10^, (c) = 8 x 10^, (d) and 
(dl) = 4 X 10^, (e) and (el) rid = 4: x 10^- Full lines corresponds to 5p^ = ±1, dashed lines to 
equilibrium initial conditions for the BCs. 
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